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Abstract. We explore possible asymmetric dark matter models using CP violating scatter¬ 
ings to generate an asymmetry. In particular, we introduce a new model, based on DM 
fields coupling to the SM Higgs and lepton doublets, a neutrino portal, and explore its UV 
completions. We study the CP violation and asymmetry formation of this model, to demon¬ 
strate that it is capable of producing the correct abundance of dark matter and the observed 
matter-antimatter asymmetry. Crucial to achieving this is the introduction of interactions 
which violate CP with a dependence. 
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1 Introduction 

Measurements of the CMB and BBN show that the baryon-to-entropy density ratio in our 
universe [1, 2], 

Yb = (0.86 ±0.01) X 10“^°, (1.1) 

is much higher than the SM predicts [3, 4]. Further, observations show us that the Standard 
Model (SM) only describes 15% of the matter in the universe, the remainder being dark mat¬ 
ter (DM). In light of this, it is very suggestive that the mass densities of dark and visible 
matter are so similar, with ^dm — where Dbm {^ vm ) is the DM (visible matter) 

density expressed as a ratio to the critical density. In many standard DM scenarios this must 
be taken as a coincidence. On the other hand, in asymmetric dark matter (ADM) models 
the DM density is due to an asymmetry which is linked to the baryon asymmetry in some 
way [5-7]. 
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To do this requires that the DM and visible sectors be related - possibly at very high 
energies. Suppose there is a conserved global symmetry between dark and visible matter such 
that B — L — D is conserved, where B and L are the usual baryon and lepton numbers and 
D is the dark matter charge. For asymmetries to form, B — L + D must be broken. Then if 
is due to a D asymmetry, the DM mass must lie in a narrow range to ensure the correct 
abundance [6], 

moM — Qd X (1.7 - 5) GeV, (1.2) 

where Qd is the charge of dark matter under D. The value of 1.7 GeV holds for totally asym¬ 
metric dark matter with Qjj = 1 where baryogenesis occurs entirely before the electroweak 
phase transition (EWPT), as will be considered in this work. 

As in normal baryogenesis, a process must satisfy the Sakharov conditions (particle num¬ 
ber violation, G and GP violation and thermal non-equilibrium), to be capable, in principle, 
of forming a particle number asymmetry [8]. Decays of heavy particles, first order phase 
transitions and the Afheck-Dine mechanism have all been proposed as mechanisms to create 
a baryon asymmetry and related DM asymmetry. A full list of references can be found in 
reviews [5-7]. In this paper we will instead study the creation of an asymmetry through GP 
violating scatterings. 

The possibility of using GP violating scatterings or coannihilations to generate ADM was 
first proposed in [9], the idea of using asymmetric freeze-out was introduced in [10] and baryo¬ 
genesis via annihilating particles was also studied in [11]. The details of such a mechanism in 
ADM were further explored in a toy model context in [12]. A baryogenesis mechanism, using 
the neutron portal effective operator, in which GP violating scatterings typically dominate 
over decays, was discussed by the present authors in [13]. GP violating scatterings have also 
been studied in leptogenesis, in which they are negligible compared to decays outside of the 
weak washout regime [14-16]. The aim of the current paper is to explore for the first time 
viable and UV complete ADM models which use the mechanism of GP violating scatterings 
during freeze-out. The DM will be stabilised using a Z 2 symmetry rather than couplings with 
a strong temperature dependence as in [10]. 

GP violating 2 0 2 processes are also crucial to the WIMPy baryogenesis scenarios, in 
which the baryon asymmetry is generated by the annihilations of weakly interacting massive 
particles (WIMPs) during freeze-out. However, the DM density is not due to an asymmetry 
in those scenarios [17-21]. GP violating scattering has also been studied in an ADM context, 
albeit for freeze-in scenarios [9, 22, 23]. 

The structure of this paper is as follows. In Section 2 we introduce the neutrino portal. 
We couple LH to an exotic scalar and fermion, the lightest of which will form DM, and 
discuss possible UV completions. In Section 3 we analyse asymmetry production in the 
minimal inert two-Higgs-doublet UV completion of the neutrino portal. We find it cannot 
produce an asymmetry of the observed size. Due to this we discuss a simple extension in 
Section 4, introducing an additional Higgs doublet, and show the observed asymmetry can be 
obtained due to the additional GP violation present. We then conclude. Alternative operators 
are discussed in Appendices D and E. 
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2 The Neutrino Portal 


When trying to write down a model of ADM, it is simplest to start with the lowest mass 
dimension gauge invariant operator involving lepton or baryon number, LH {H r\j (1,2,-!)). 
This must be paired with at least two non-SM particles in order to stabilise dark matter 
with a symmetry, leading to a neutrino portal [24, 25]. This can be achieved via the effective 
operator 

gLYHcP, (2.1) 

where gf is a coupling with inverse dimension of mass, L (H) is the SM Lepton (Higgs) doublet 
and Y ((/>) is an exotic fermion (scalar). The latter two fields will comprise our dark sector, 
in which the lightest particle’s stability will be ensured by a I 12 symmetry. Both cp and Y are 
odd under this symmetry, and all SM particles are even. 

It is preferable to have only one dark sector particle carry a baryon number equivalent 
D, because if both particles did, then they must both be stable. We can define a dark baryon 
number D = {N^j, + Ny)/‘2 (where Nx is the particle number of X), so that B — L — D is 
conserved and B — L + D is broken, as desired for ADM. Unfortunately if both cj) and Y carry 
D we have a U(l) symmetry in addition to B — L — D conservation. This is obtained by 
rephasing only the dark sector particles, giving 


A(/)-Ay = 0, (2.2) 

where XX = Nx — Nj^. If (/> or T decays, (2.2) requires that the asymmetry stored in the 
remaining particle vanishes.^ Then baryogenesis does not occur as B — L = 0.^ This problem 
is avoided if only one dark sector particle carries D, e.g. D = 

To keep both (j) and Y stable would require kinematically disallowing decays with a rather 
implausible mass difference smaller than the neutrino masses. This mass difference cannot be 
obtained as one particle species must go out-of-equilibruim above the EWPT (so sphalerons 
can reprocess AL into AB) and DM must be of order a few GeV. While it is possible that 
additional interactions could prevent Y 01 (j) from carrying an extra dark conserved charge, 
for simplicity we will not consider that scenario. This leaves two possibilities: either (j) is a 
complex scalar and Y a Majorana fermion, or (f) is a real scalar and Y a Dirac fermion. We 
choose the former case as it leads to more CP violating phases and permits a simpler method 
to annihilate the symmetric component. Asymmetry production should be the same in both 
cases, though they will have phenomenologically different dark matter (scalar versus fermion). 

We envision the Ts as the most massive particles (barring the mediator), going out-of- 
equilibrium before our dark matter candidate (/>. As they do not store an asymmetry, they 
can safely decay without removing the asymmetry generated by the scatterings. To satisfy 
unitarity constraints, which we will discuss in Section 3.1, there must be two generations of 
Ys and one of cj). Our ADM will be (j)-, which carries a non-zero baryon number equivalent 
D = N^, preserving B — L — D. 


there are additional dark sector particles carrying D, then an extension of this argument still holds: at 
least two must be stable. As this is non-minimal, we do not consider that scenario in this work 

^ As a special case this can be avoided if the scatterings freeze-out before the EWPT but the particle decays 
after the EWPT, because B and L are no longer related by electroweak sphalerons. 
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Figure 1. The interaction LYHcj) together with its UV completions. We label them as case 1 (top 
right), case 2 (bottom left) and case 3 (bottom right). 

In this paper we consider UV completions of the effective interaction, eq. (2.1). The 
simplest tree level possibilities are depicted diagrammatically in fig. 1. The thermal histories 
of all UV completions are essentially the same. At high temperature the Vs are kept in 
equilibrium by the rapid 2 —)• 2 interactions. These interactions eventually freeze-out and 
the Vs depart from thermal equilibrium, satisfying the third Sakharov condition. As these 
scatterings violate CP, a non-zero B — L is generated along with an associated D asymmetry 
stored in the 4>, creating a baryon asymmetry via the sphalerons (as long as this occurs above 
the EWPT). Subsequently, the Vs decay into leptons and cj), potentially leading to further 
asymmetry formation. To see which process dominates, as well as the size of the asymmetry 
formed, it is necessary to solve the Boltzmann equations, which we will do in Section 3. The 
three simple UV completions are as follows. 

• The first case requires the mediator, /, to be a vector-like lepton. The Z 2 symmetry 
prevents / from Yukawa coupling via the Higgs to the SM leptons, avoiding most of 
the limits on vector-like leptons [26, 27]. For CP violation to arise at one-loop it is 
necessary to include two copies of /. 

• The second case has a Dirac SM singlet, / as intermediary. The Z 2 symmetry precludes 
the V from acting as an intermediate particle and coupling to the leptons, which would 
destroy the global B — L — D. If V is Dirac, rather than Majorana, this operator is 
very similar to [9]. In [9], a Majorana fermion mediates an interaction between LH and 
a mirror copy, L'H' , using different initial temperatures to create out-of-equilibrium 
conditions. CP violation is the same as in the first completion. 

• The third case is an extension of inert two-Higgs-doublet models (IDM) [28]. The 
intermediate particle, H 2 , has the same quantum numbers as the SM Higgs but cannot 
play the same role. While CP violation is possible with just one inert Higgs doublet, 
we will show that to get a sufficient asymmetry to form, two inert Higgs doublets are 
necessary. While regular IDMs do have a dark matter candidate, they do not explain 
the ratio of dark and visible matter; the mass ranges required for H 2 to be ADM instead 
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of (j) are excluded by a combination of collider searches, electroweak precision tests and 
direct detection constraints [29, 30].^ 

As the IDM completion is particularly illustrative we will study this case in more detail. 

2.1 Inert Higgs Doublet Completion 

To see a specific model in action, we will study the IDM completion of the neutrino portal. 
As in IDM, we add a massive scalar SU{2)l doublet which will function as our heavy inter¬ 
mediary. We also introduce a new parity, with the Lagrangian symmetric under H 2 —)• —H 2 , 
Y —)• —Y, and (p —)• —cp. All SM particles are even under this parity. We call H 2 inert as it 
does not acquire a non-zero vacuum expectation value (VEV), and does not play a role in the 
fermion mass generation. Where this work differs from traditional IDMs is that in our case 
H 2 carries a non-zero B — L.^ The relevant additions to the SM Lagrangian are 

AC = -XiaH2LiYa-KHiHl^-Xi\Hi\^\H2\^-X2\HlH2\^-Xs\Hi\^\cP\^-X4\cP\‘^ + H.c. (2.3) 

This induces (2.1) after integrating out the mediating H 2 . After considering field rephasings, 
there are three physical CP violating phases, that we write as 


'i = Arg(Ai 2 ), 

(2.4) 

>2 = Arg(A 22 ), 

(2.5) 

I 3 = Arg(A32), 

( 2 . 6 ) 


with all other couplings real (without loss of generality). For our example solutions we choose 
= 02 = O3 = 7 r / 4 . 

We consider the Y to have a mass higher than the scale of the EWPT, and H 2 to be 
heavy enough to be considered always off shell, typically several orders of magnitude heavier 
than the next lightest particle, Y 2 . As the neutrino portal does not violate baryon number 
directly, electroweak sphalerons are required to reprocess the lepton asymmetry into a baryon 
asymmetry. For the parameters considered in this paper, all relevant processes finish before 
the EWPT. 

We assume that there is an additional process that keeps the (p in equilibrium and anni¬ 
hilates the symmetric component, but remain agnostic as to the details of this process. For an 
example, cp quartic coupling to a light real scalar field could be used for this purpose. In more 
baroque extensions it is also possible for cp to decay into lighter particles, and for them to have 
new interactions that annihilate their symmetric component. While for the purposes of this 
paper the exact method is irrelevant, different mechanisms will provide different detection 
prospects. For example, if (p annihilates into a light, stable state then there is a significant 
contribution to dark radiation [ 6 ]. 


^An exception to this can be found in [30], where an inert Higgs transfers a pre-existing asymmetry between 
the dark and visible sectors. 

'^The main effect of this is on neutrino masses. In standard IDMs neutrinos gain Majorana masses through 
radiative corrections [31], which requires a term As (HI Hi )^. As this term breaks B — L — D, for simplicity we 
consider instead neutrinos gaining a mass through the Higgs mechanism (with the addition of right handed 
neutrinos even under the Z 2 ). This can be a Dirac mass or a Majorana mass if lepton number is softly broken. 
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Originally it was thought possible to use CP violating scatterings to simultaneously 
create an asymmetry and annihilate the symmetric component [10]. Unfortunately, this is 
impossible as it would require the freeze-out temperature to be of order mDM/25. For op¬ 
erators like those of [10], decays prevent an asymmetry forming at arbitrarily low freeze-out 
temperatures. If we make the coupling strength large enough to delay scattering freeze-out 
until Tfreeze-out ~ ?^dm/ 25, the lifetime becomes short enough that decays and inverse decays 
keep the particle in equilibruim, so Sakharov’s conditions are never met. We hnd numeri¬ 
cally that the lifetime becomes shorter than tfreeze-out well before Tfreeze-out ~ w,dm/25. As 
discussed in Appendix D.2, models that might have this feature are disfavoured by limits on 
non-SM coloured particles. 


2.2 Phenomenological Constraints 

Due to the mass range of H 2 we are considering most of the couplings related to asym¬ 
metry formation in our theory are unconstrained, but there are still some restrictions from 
electroweak considerations.^ 

2.2.1 Electroweak Phase Transition 

While the population of H 2 is negligible during the EWPT, it is still worth considering the 
effects of the EWPT on H 2 as it couples to cj). In addition, while <p is a SM singlet it can 
still potentially influence the EWPT. During electroweak symmetry breaking (EWSB), Hi 
acquires a VEV v and acts as the SM Higgs; on the other hand H 2 only experiences mass 
splitting. It is possible to parameterize H 2 as 


H2 = 


-”2 


(2.7) 


where and H are complex scalars. After EWSB we have 

m?^- = -h Xiv^. 


( 2 . 8 ) 

(2.9) 


There is also a contribution to the mass of </> when Hi acquires a VEV, as well as mass mixing 
between H 2 and cj). This gives us a mass matrix in the ((/>, H 2 ) basis 


KV 



( 2 . 10 ) 


We expect this to diagonalise to a heavy state and a light state. We will label the light state 
(f) as the dark matter admixture is mostly (f), with mixing angle Ku/m^o- We have a light 
enough dark matter candidate for small mixing angles, 

vkIitiho < {ml + Aan^)^/^. (2.11) 

®While directly searching the full possible mass range of H 2 is impossible at the LHC, it is possible to rule 
out high scale baryogenesis in general. If low scale lepton number violation (LNV) is observed, either through 
direct observation of LNV at the LHC or through a combination of neutrinoless double beta decay and lepton 
flavour violation, then the washout induced by this LNV makes any high scale baryogenesis irrelevant [32, 33]. 
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This requirement for small mixing is the only real constraint relevant to asymmetry forma¬ 
tion. For K ~ my and ~ 0, we have m^, ~ X^v^. As we envision cj) to be ADM, we must 

have As < 5 X 10“^ to get the correct relic density ~ 1.7 GeV as in (1.2)).® 


Introducing (p can potentially alter the EWPT, but only via its quartic coupling to 
Hi. One might worry that if the EWPT is made first order, there would be two competing 
methods for baryogenesis - electroweak baryogenesis and scatterings. From the analysis in 
[34] it can be shown that the EWPT is first order if 


^ - -^Xjj.^iXs + A 4 ) > 


1287r2 


Ah, I ^ + 1 ) - ^ 


m 


( 2 . 12 ) 


where Xh^ is the quartic self coupling of Hi and yt is the Yukawa coupling of Hi to the top 
quark. As A 3 is small, the EWPT remains second order. We make the approximation that 
the electroweak phase transition occurs the same way as in the SM. 


2.2.2 Electroweak Precision Tests 

As H 2 still couples to the electroweak gauge bosons, electroweak precision tests are, in prin¬ 
ciple, sensitive to H 2 . More precisely the electroweak precision tests are only sensitive to the 
mass splitting of H 2 [28]. Standard electroweak precisions tests can be cast in terms of the 
"oblique" or Peskin-Takeuchi parameters S, T and U which parametrize the radiative correc¬ 
tions due to new physics [35]. From [28] the main contributions to electroweak precision tests 
(at one loop order) can be written as 


AT ~ 

A5~ 


A2 


327r2Q; 

X2v‘^ 




0.4A2, 

~ 8 X 10"® 


X A2 


10 TeV 


m. 


h: 


(2.13) 

(2.14) 


The experimental values are ATjj^q = 0.08 ± 0.07 and A5{/=o = 0-05 ± 0.09, so A 2 < 0.2 is 
required to satisfy electroweak precision tests [36] . Fortunately, A 2 does not affect asymmetry 
production, so this constraint does not seriously affect the model. AS is negligible for all 
sensible parameter choices. 


2.2.3 Self Interactions 

While the self couplings of p do not affect asymmetry production directly, they will contribute 
to the thermal mass of p. The most stringent limits on dark matter self interactions come from 
the Bullet Cluster and similar colliding clusters. From this there is the constraint, [37-39] 


Al < 4 X 10® 




1 GeV 


(2.15) 


This leaves A 4 essentially unconstrained. As the thermal mass of p will be determined by A 4 , 
this freedom is useful for the kinematics. We use the thermal masses described in Appendix B 

®In this case A 3 is much too small for (j> to be detected at the LHC via invisible Higgs decays. However, if 
there are significant cancelations between mp and it is possible for A 3 to be large enough for the invisible 
decays of the SM Higgs into 4> to be detectable at the LHC, though this is not required for our scenario to 
work [40]. 
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throughout this paper. 


Fortuitously only one of these constraints affect parameters required for asymmetry 
formation, so we are free to choose parameters to maximise the asymmetry formed. The only 
real constraint is (2.11). 


3 Finding the Asymmetry 


3.1 Notation and Tools 

Now that we have invented some interactions satisfying the Sakharov conditions, we need to 
know the magnitude and nature of the asymmetry they generate. Our main technique for 
calculating the evolution of an asymmetry will be Boltzmann equations. We adopt the same 
W (V’, a ...—)• i,j ...) = n^ria ^a ...—)■ ...)) notation as [13], where a ... — 

i,j ...)) denotes a thermally averaged cross section, and parameterize the CP violation as 


^ iT(V^,a... i, j ...) - VF(i, j ... -» V;,a...) 

a... ^i,j ...) + W{i,j Ip, a...)' 

We define the CP symmetric reaction rate density as: 


fPsym — 


W{i,j ^ pj, a...) + W{'ip, a.^ i,j ...) 


(3.1) 


In general both Wsym will be temperature dependent. As we will be using 

Maxwell-Boltzmann statistics as a (good) approximation throughout this paper, we can factor 
out the chemical potential from the phase space 




(3.2) 


where refers to the equilibrium values when the chemical potential is zero, and define: 


_ 




(3.3) 


Putting this all together, for the process ip + a + b . —> i + j + ..., evolving in the 
Friedmann-Robertson-Walker metric, the evolution of is given by [41]: 


— -F ZHn = Ciip), 

where H is the Hubble expansion rate and the collision term C{'ip) is: 


(3.4) 


C{lp) — ^Vsym 


TiTj ... (1 £i/i,a...—>1,^...) •••(IT £'!/i,a...—j...) 


(3.5) 


We compute {va{'ip T a —t i T j)) from the cross sections by using the result from [42], 

{vaipj + a^i + j)) = J p^aE^I,EaV,el(TKi ds, (3.6) 

where s is the square of the centre-of-mass energy, pr^jja is the centre-of-mass momentum of 
Ip and a and Ki{x) is a modified Bessel function of the second kind of order i. As a note 
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of caution, we will not use the Einstein summation convention when discussing Boltzmann 
equations. 

The requirement that CPT and unitarity must hold imposes an important restriction 
on our CP violating terms. In our language the unitarity constraint is given by: [43-45] 

^ W{i ^j) = Y, W{j = = Wii ^ j), (3.7) 

j j j j 

where i is the CP transform of i. In general (3.6) must be solved numerically, which can allow 
small artificial violations of unitarity to arise. To deal with these small errors, we enforce the 
unitarity relations (3.7). 

To calculate the CP violation we use Cutkosky rules [46]. It should be noted that these 
are the cutting rules for zero temperature quantum field theory; when thermal effects dom¬ 
inate it is necessary to use the rules contained in [47, 48]. As we will only be interested in 
the dominant contributions to CP violation near freeze-out, we will use the zero temperature 
cutting rules. 

With these techniques, it is possible to analyse the asymmetry production. 

3.2 Interactions and CP violation 
3.2.1 Scatterings 

We now catalogue the relevant interactions. Our CP violating scatterings are given by: 

W{YaL ^ Fk/)) ^ LYa) = (1 + ea)Wa, (3.8) 

W(LYa ^ ^ YaL) = (1 - ejWa, (3.9) 

where we have included an implicit sum over the three families of leptons. Because all the 
leptons have essentially the same mass and chemical potential, it is not really necessary to 
consider them as separate species except when summing the couplings for a process. These 
annihilations will be the main generators of the asymmetry. The relevant unitarity constraint, 
derived from (3.7), is 

eilTi = — 62 ^ 2 - (3.10) 

The only CP violation for the process YaL Hi4> comes from interference between 
graphs in fig. 2, involving a Majorana mass insertion. No other one loop graphs lead to a 
complex combination of couplings. Taking generic couplings A ~ Xia for the Lagrangian in 
(2.3), the tree level cross section scales as 


av ~ 




(3.11) 


Denoting a as the CP conjugate cross section, the CP violation at the cross section level 
scales as 


{a — a)v 


X^K‘^mYi'mY2 


(3.12) 
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Figure 2. Graphs contributing to Ca- The one loop graph (right) has a Majorana mass insertion. 


Inserting this into (3.6), noting the momenta and energies scale as ^/s and the Bessel function, 
which falls as Ki(z) ~ ^ 1; imposes an effective cut-off on the integral at 

^/s ~ T, one finds 


IFi oc 




(3.13) 


and 


eiVFi oc 


^4.2 




'IpG 


(3.14) 


Hence, from this dimensional analysis, it is clear there is no T dependence in the CP violation, 
ei for T > My2-‘ We also have CP conserving scatterings, which we label 


W{YaHi ^ Let)*) = Ta, 

(3.15) 

W{YaCt> ^ HlLi) = Ua, 

(3.16) 

W{YaL ^ Y,L) = Sab, 

(3.17) 

W{LL ^ YaYb) = Pab- 

(3.18) 


While LL —>■ YaYf, is not technically CP conserving, CP violation in this term only leads to a 
flavour asymmetry. As Y is Majorana, and not ultra-relativistic, there is no method to store 
this asymmetry or transfer it to a non-zero B — L so these terms can be safely neglected. For 
the t-channel graphs, no loop graph with an absorptive part is kinematically allowed and so 
there is no contribution to the CP violation. The expressions for the cross sections of the 
above processes are given in Appendix A. 


3.2.2 Decays 

We consider the following decays of Yf, I 2 and Hi, though we are only interested in the CP 
violation of the Y decays. We will label these decays by: 


T{Hi ^ YaLcj)*) = 

(3.19) 

r(y2 ^ YiLL) = T2A 

(3.20) 

r{Y2^LHiej)) = i(l + e^)r2s 

(3.21) 

r{Y2 ^ LH*i(l)*) = l(l-e^)r2s 

(3.22) 

r(yi ^ LHi^) = r(yi ^ LH*i(j)*) = ^Pi, 

(3.23) 


^This is, of course, ignoring the implicit temperature dependence of the thermal masses, but the overall 
point still stands. 
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Figure 3. Graphs contributing to ea- The one 





graph (right) has a Majorana mass insertion. 


where we have parameterized the CP violation in decays by eo- Note that (3.19) is kinemati¬ 
cally allowed at high temperatnres dne to the thermal mass of Hi. Contribntions to come 
from the graphs in hg. 3. We can see from this hgnre that Yi decays are CP conserving: CP 
violation wonld reqnire an on shell Y 2 in the loop, which is kinematically forbidden. There 
is a farther contribntion to the CP violation from the 3 —)> 3 scattering LHicf) —)> 
with the real intermediate Ys snbtracted [14]. While the CP symmetric component of this 
scattering is negligible, CP violation in this process serves to cancel the CP violation in the 
decays at thermal eqnilibrinm [13]. This can be seen by applying the nnitarity relation (3.7) 
to the state LHicj), to derive 






(3.24) 


which we inclnde in onr Boltzmann eqnations. 


3.3 Boltzmann Equations 

The chemical potentials of the SM helds and (p depend only on the B — L asymmetry (see 
Appendix C), so we have only three conpled differential eqnations to solve. These Boltzmann 
eqnations are similar to those in [13]; in both cases we have two heavy Majorana particles 
which experience CP violation throngh decays and 2 —)• 2 scatterings. For Yi and 12) we 
have: 


dUYi 

dt 


+ 3HnYi = Tin^_^ r-fu^r^)/2 - rY^ - r2Any^ 


rirjTY., - rY2 




+ Wi 


- eiWi 


rHH4> - (rj - n) 


+ Ui 
+ 2 Pii 


- VY, (r 0 + r^) 


+ Ti 
— *S'i 2 


rBir^ + - ty^ {rj + n) 

irj + ri){rY^ - ryj 


rjri - rY^ 


+ -P12 


rjn - rY^VY^ 


(3.25) 
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and 


dnY2 

dt 


+ 3HnY2 = {nrjj^r-^ +rjrH^r^)/2 - ty^ + 




+rH^- nr^'rY2 - rjr^rY2 


rirfY^ - rY2 
eq 


+ W 2 
+ T 2 
+ 5'i2 


rH^r^ + - ry^ (rj + n) 

- ^y2 


+ eiT^i 


- ^€Dr2BnY2 - rjrn.r^ 

rn^r^ - rjj^r^ + ry ^(rj - n) 


+ U 2 


^H^rj + rjj^ri - ry^(r^ + r^) 


(rT + n)(ryi -ryj 


+ 2 P 22 


rjn - rY2 


+ P 12 


rjri - ry.ryj 


(3.26) 


The Boltzmann equation for B — L is: 
dns-L 


dt 


+ 3HnB-L 


= Tin^y 


- rjrH^r^ 


+ T2Bn^. 


nrjTj^ - rjVH^r^ 


+ r^H'j-s.yjn®'^ 


Hi 

eq 


"^Hi - - rir^VYi + rjr^rY^ 

+ THi^Y2nH^ rui -rH^ + rtr(j,rY2 - rir^Y2 + rjr^rY2 


- 7,^D'^2BrfY2 


2rY2 - {rirjj^r^ + rfHir^,) 


+ eiWi ^{ri + r;)(ry2 - ryj 
+ Y1 ~ ~ ^HiTYa - rn^Ya + r^j - 


a=l,2 


+ Ua 


- rniTj + r^VY, - r^VY, ) 


driD 

dt 


+ SHtid- 


(3.27) 


The terms with ei and eo are source terms for the asymmetry, while all other terms tend 
to wash out the asymmetry. We have used (3.10) to write the Boltzmann equations solely 
in terms of these two CP violating terms. The decay rates which appear here are thermally 
averaged (see Appendix A). 


We can now numerically solve these to see the evolution of an asymmetry. 

3.4 Numerical Solutions 

Our numerical solutions to the Boltzmann equations are stable under changes to the preci¬ 
sion, starting temperature (for initial temperatures T ^ HIY 2 ) and initial conditions. 

We start our solutions at high temperature, where 2 —>■ 2 scatterings can wash out any pre¬ 
existing asymmetry, and then track the evolution of i? — L down to the EWPT. To good 
approximation, at the EWPT sphalerons simply switch off, freezing the value of B. In gen¬ 
eral, if Yf or I 2 is sufficiently long lived they may come to dominate the energy density at 
some time, causing the universe to become matter dominated (instead of radiation dominated 
as we assume). This would lead to a dilution factor due to reheating when the Y decay [49]. 
For the interesting regions of parameter space (those which lead to large asymmetries) this 
condition is not satisfied but for completeness we do include the dilution factor in our code. 
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Figure 4. Asymmetry formed as a function of k and m/Zj. There is a ridge of values where the 
asymmetry formed is significant, corresponding to a freeze-out temperature of order My^/S. This 
temperature is due to the interplay between decays and scatterings: increasing the couplings leads to 
'''Y ^ tfreeze-out, which forccs the Y to remain in equilibrium, and decreasing the couplings decreases 
CP violation. Maximal asymmetry corresponds to At ~ and to/Jj ~ 10^ x my^. In this example, 
myj = 90 TeV and my^ = 100 TeV. 


For regimes where the two Y are relatively close in mass, scatterings dominate the asym¬ 
metry, as also found in [13]. As there are two mass scales in this problem (once we choose 
a mass for Yi and Y 2 ), we scan over k and keeping in mind (2.11). The asymmetry is 

maximised when k ~ and mjyj ~ 10^ ™Yi ■ Unfortunately, this model does not generate 
the full asymmetry required for baryogenesis; from fig. 4 it can be seen that the maximum 
asymmetry is of order Yb ~ 10“^^. 


This is quite puzzling: with such similar Boltzmann equations and asymmetry produc¬ 
tion methods, how can this model fail where [13] succeeded in obtaining Yb rsj 10 The an¬ 
swer lies in the CP violation. Whereas in the neutron portal case there were graphs that gave 
the CP violation in scatterings the property e{Xu —)• dd)W{Xu —)> dd) oc T‘^W{Xu —>■ dd), 
there is no temperature dependence in our CP violating terms. As was suggested in [13], 
which we will show explicitly, in the neutron portal case this temperature dependence made 
the scatterings relevant at a higher temperature, enhancing the asymmetry production. To 
obtain similar temperature dependence, we must have CP violating graphs without Majorana 
mass insertions. To this end we must include a second copy of the mediating scalar. In fact, 
this conclusion also holds for [12, 13]: to get the full asymmetry found in the EFTs studied it 
is necessary to have two heavy intermediate scalars regardless of the number of CP violating 
phases. 
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Figure 5. New one loop graphs contributing to CP violation. The graph on the left (right) contributes 
to €a (cd)- Neither graph has a Majorana mass insertion. 


4 Extended Neutrino Portal 

A simple extension of the neutrino portal, which allows the generation of the observed Yg, 
is the addition of another inert Higgs. Three-Higgs models have been explored in both an 
inert and general context [50, 51]. We will not write down the full potential, as most of the 
terms are irrelevant for asymmetry formation, but simply note that to satisfy electroweak 
precision tests it is necessary to avoid significant mixing between the two inert scalars [52]. 
The Lagrangian is now 

A£ = -mjjjHpf - XiapHpLiYa - KpHiHj^cf) + H.c, (4.1) 

where p = 2,3. With 14 relevant couplings there are now nine CP violating phases which, for 
the sake of the example, we will choose to be 

Arg(Ai2i) = 0 , Arg(A22i) = Arg(A 32 i) = Arg(Ai22) = Arg(A2i2) = ^ 

Onr nr ^TT TT 

Arg(A222) =—, Arg(A3i2) = —, Arg(A322) = —, Arg(K2) = -• (4.2) 

5 10 5 4 

We have chosen phases that avoid cancellations between the various CP violating graphs. 
There are now additional graphs which lead to complex couplings, involving a closed fermion 
loop (fig. 5). As there are no Majorana mass insertions, these graphs exhibit the desired 
temperature dependence, eiWi oc T^VPi (see fig. 6). This can be seen by looking at the 
second term in (A. 12), which corresponds to this new graph; relative to the first term, which 
is due to the original CP violation with Majorana mass insertions (fig. 2), there are two 
extra factors of energy and momentum, E and p. When velocity averaged with (3.6), at 
high temperatures E ^ p ^ T. At low temperatures the CP violation becomes constant, 
though due to the different combination of phases involved it is not necessarily equal to the 
CP violation in fig. 2. There are similar contributions to the CP violation in decays, which we 
also include. Armed with this new CP violation, we can again solve the Boltzmann equations 
for the evolution of B — L. 

4.1 Numerical Solutions 

With the temperature dependent CP violation, we see a significant increase in the asymmetry. 
As long as there are no significant cancellations between the various contributions to the 
the full baryon asymmetry of the universe can be generated (see fig. 7). Fittingly for ADM, 
cancellations are minimised when the Xiap are not symmetric: it is preferable for the leptons 
to couple more strongly to one of the generations of the Y and for not all generations of 
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Temperature [GeV] 


Figure 6. CP violation in the extended neutrino portal with my^ = 90 TeV and = 100 TeV. 
At high temperature the CP violation from the graphs in fig. 5 grows as T^, and the CP violation 
from the graphs in fig. 2 is constant. At low temperatures both contributions become constant, 
and approximately equal in value. We have plotted both sources of CP violation with phases and 
couplings chosen to maximise their CP violation, though no parameters maximise both simultaneously. 
At T ^ 3 X GeV, where maximal asymmetry production occurs, the CP violation from the graphs 
in fig. 5 dominates. 


leptons to couple with the same strength. This asymmetry in the couplings does not need 
to be more than an order of magnitude for significant enhancement, as in fig. 8 . Similarly, it 
is preferable for K 2 and K 3 to be an order of magnitude different, and mjyj and to be 
within an order of magnitude of each other. Comparing this to the asymmetry formed when 
only the mass insertion graphs are included (fig. 7) we see that the asymmetry starts forming 
earlier, culminating in a significantly higher asymmetry. 

Due to the similarity of the asymmetry production, Boltzmann equations and CP viola¬ 
tion between this model and [13], it is clear that main asymmetry production in the neutron 
portal EFT studied in [13] was also due to this temperature dependent CP violation (rather 
than the mass insertion diagrams which are also there in the neutron portal case). This pro¬ 
vides compelling evidence that for scatterings to dominate over decays a heavy intermediate 
particle is necessary (to provide the dimensionality for temperature dependence). As many 
cosmological models are EFTs, this will often be satisfied. Further, multiple intermediate 
particles (leading to bubble graphs) seem to be a generic feature of scattering models, being 
necessary in not only this work but also [12, 13, 17]. 


5 Conclusion 

We have shown that it is possible to create a realistic UV complete model of ADM that uses 
scatterings to generate the asymmetry. Further, we have demonstrated explicitly the impor¬ 
tance of temperature dependence in the CP violating terms. 

To do this, we introduced a new model, using an neutrino portal. By examining the UV 
completions, we separately studied temperature independent and temperature dependent CP 
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Temperature [GeV] 


- Y^y, 

- Yy, 

- FV, 

- Yy^ 

l^fileocr" 

\YB\eocT^ 

I Decays 


Figure 7. Example solution with wyi = 90 TeV, my^ = 100 TeV. The asymmetry generated is 
Yb = 1.0 X 10“^°. We use e oc T° and e oc to refer to the CP violation in the scatterings from 
figures 2 and 5, respectively. The asymmetry generated by the decays is negligible. 




2x10^ 


4 X 10^ 6 X 10^ 8 X 10^ 

Mass of H 2 [GeV] 


1x10*^ 


Figure 8. Asymmetry formed vs niH^i with mu^ = 2.bmH2- asymmetric couplings the asym¬ 
metry formed is as much as two orders of magnitude higher than when all the couplings are equal. 
Asymmetric couplings can avoid cancellations between the various CP violating graphs. The asym¬ 
metric couplings have been chosen for this example to maximise asymmetry production. Coupling Yi 
and Y 2 to one generation of leptons, with the couplings differing by an order of magnitude, is sufficient 
for a large asymmetry to form. When there is little or no discrimination between the couplings, the 
asymmetry formed is essentially the same as the unextended neutrino portal. If the two inert Higgs 
have the same masses and couplings then the temperature dependent CP violation cancels exactly. 


violation, with the latter proving to be far more signihcant. For fntnre ADM model bnilders 
to see a signihcant asymmetry cansed by 2 —)• 2 scatterings, qnadratic temperatnre depen¬ 
dence in the CP violation will be a key featnre. 
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There are many other potential operators if one is willing to allow a slightly more compli¬ 
cated dark sector, inclnding variations on the WIMPy baryogenesis operators. An interesting 
avenne for fntnre research wonld be the annihilation of the symmetric component of the (j) 
nnmber density; while the same scatterings cannot simnltaneonsly create an asymmetry and 
annihilate the symmetric component of ADM in existing models, other minimal options can 
still be explored. 
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Appendix 

A Cross Sections of the Neutrino Portal 

For the sake of completeness, here we catalogne the cross sections of the case stndy in Sections 
2.1 and 4. Ei, E 2 , E^ and FI 4 refer to the energies of the particles in the order listed. For the 
2 —>• 2 scatterings, we will denote the initial momentnm pi and the final momentnm pf. All 
cross sections are written in the centre of mass frame, except when otherwise stated. 

A.l Minimal Neutrino Portal 

A.1.1 Cross sections (CP conserving component) 

We inclnde some of the factors from (3.6) to highlight the symmetry nnder interchange of 
particles. 

YaLi ^ Hict> 



1 2^2 1 

PiEiE 2 av= “4 piPf{EiE 2 +pi). 

iOTT Y 

(A.l) 

E 4 > ^ YaHi 

1^2x2 1 

PiEiE 2 av = piPfEiEs. 

STTy/smjj^ 

(A.2) 

L^Hi ^ Ya(^ 

1^2x2 1 

PiEiE 2 av = piPfEiEs. 

STTy/smJj^ 

(A.3) 

YaLi ^ YfyLj 



PiEiE 2 av = 

|A^ 1 

, PiPf [{E1E2 +pj){EsE^+p}) + EiE 2 E^E, + l/ 2 p 2 p 2 ] , 

47 r srrijj^ 

(A.4) 

LiLj YaYf, 




I \ 2 \2 I I I \ 2 \2 I 

PiEiE2av = V PiPf [E1E2E3E4 + l/ 2 p 2 p 2 ] . (A. 5 ) 
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A.1.2 Cross sections (CP violating component) 

Y2Li ^ Hicl, 


PiEiE 2 {a - a)v = 


mY^rnY2 Im(AjiAjiA*2A*2)«" 


3 


Glvr^ sm 


6 

H2 


■PiPfPloopE2EL^, 


(A. 6 ) 


where Ej is the energy of Lj in the loop and pioop is the momentum of the (on shell) particles 
in the loop. The factor of myjmy 2 comes from projecting out the Majorana masses. CP 
violation in YaE —?■ Yi,Lj is similar but only leads to flavour violation and is so neglected. 


A.1.3 Decays (CP conserving component) 

To calculate the three body decays, we used the usual trick of decomposing N-body phase 
space into a series of 2-body phase spaces. All integrals are numerically integrated. The decay 
rates appearing in section 3.3 are in fact thermally averaged [14], 


-pthermal Ai (m/T) 

“ K2{m/T) 


(A.7) 


Since the relative masses of the particles changes with temperature, at different times the Y, 
Hi and (p can all decay depending on the thermal masses (in particular the choice of A 4 ). We 
will only write down the Y decays as the others are similar and not particularly important. In 
our example solutions we chose A 4 so that at high temperatures Hi —)• LiYap occurs. Keeping 
this in mind we now catalogue the decays. 


Ya 


LiHlP* 


r= J ml [\ls + pI{s) + sj ml + pl{s^ pl{s)p^{s), (A.8) 

where El{s) is the energy of Li, pl{s) [^^(■s)] is the momentum of Lj \<p] in the centre-of- 
momentum frame [rest frame of the mediating particle 7 / 2 ] and s is p^PHifi- When we say 
the rest frame of H 2 , we mean a hctitious on-shell particle with mass ^/s in place of H 2 . 


Y 2 


LjYiLj 


(myj-mi.)- 


T = 


(rrii +myj )2 


. 2 I + |A|iAf 2 l) 


\/s+P%{s) + +P%{s)) (A.9) 


xPLi{s)pL^{s) (^Elj{s)Ey^{s) +pi.(s)) 


where El.^.^{s) and PL^y^is) are the energy and momentum of Lj q), respectively, and s is 
PyPl ii- Note that Pl (s) is in the rest frame of H 2 and all others are in the rest frame of 

hI 
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A.1.4 Decays (CP violating component) 

Fa ^ LiH*,cP* 

r - r = J ds 

Low 

(A.IO) 

where s is Pi PHni and Low = Max[(m 0 + (w-Vl + nT-Lj)^]- The kinematics of Li were 

written in the centre-of-momentnm frame, and the other particles kinematics are in the rest 
frame of Lfa- 

A.2 Extended Neutrino Portal 

When a second inert Higgs donblet is added, H 3 , we get many more diagrams contribnting 
to the varions processes. For most processes, these diagrams are simply allowing the inter¬ 
mediate particle to be H 2 or H 3 , with some interference terms. Dne to the triviality of the 
extension, we do not write these down. Bnt there are important new contribntions to the CP 
violation in both the decays and scatterings. The new cross section is: 


^ my, ^ ^ / 2 2 ^ ^ ^ ^ 

2 ^ ^4 3/2 6 ^Lj[S) [niY^ - S - TTIlJ PLi{s)pL^{s)p^{s), 


Y 2 Li ^ Hicl> 


PiEiE 2 {cr - a)v 


E 

jprs 


my^my^ Im(AiirAjlsA*2pA*2rK*Ks 


±1D 11y 


-PiPfP\oopE 2 ELj 


(All) 


A, 


+ 


lm{Xjir^jls^i 2 p^i 2 rKpKs^ 


jprs 


327r2 


sm|r 771% m% 

lip liy 11 s 


PiPfPloop {E 1 E 2 Epf) [EljEy^ Epioop) > 


where p, r, s G 2, 3 and Ey-^ is the energy of Fi in the loop. The first term is essentially the 
same as (A. 6 ), and is dne to an extended analogne of hg. 2 . The second term is dne to hg. 5: 
as there are no Majorana mass insertions we have a more complicated fnnction of energy and 
momenta in place of the Majorana masses. We note that if all the masses and conplings are 
the same for H 2 and iLs the second term cancels, and CP violation is rednced to the minimal 
nentrino portal. The contribntions to the decays are given by: 


Fo 


LiHt 


(my -niL)^ 


r-r = - 


l™(AilrAjlsAj2pA •2rKp'^s) \ / 2 2 \ 

^^2^ /o^^4.3/2^2 _2 _2 ^lAs) [my^-s-rriLj 


(2T7Ys'^Em% 771% m% TTlyo 

Low Tip tlr 71 s ^ 

X PLiis)pLjs)p^{s) 

{my^-mL% 

2 lm{X%^XjisXi 2 pXi 2 rK*Ks) 


+ 


dsY, 


{2'K)^S^/‘^m%^77l%^7n%^77lY^ 

X {El^Ey^ + pfoop) PLiis)pLjs)p^{s). 


2 ELj{s)i^^s+pljs) + ^ml.+pl^{s 

(A.12) 


Similar to the CP violation in the scatterings, in the limit of degenerate conplings and masses 
the second term disappears. This completes the catalogning of the cross sections nsed to solve 
the Boltzmann eqnations in section 3.3. 
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B Thermal Masses 


While throughout this paper most thermal held theoretic effects are neglected, there is one 
that must be included. At high temperatures particles can gain an effective mass through 
interactions with the plasma. For Hi, L and (p, these thermal masses are kinematically 
significant. H 2 is too massive to appear in the plasma at the temperatures we consider and 
so does not acquire a thermal mass, nor does it contribute to thermal masses. Similarly, 
as the Y only interact via H 2 they also do not acquire a thermal mass. The main effect 
of thermal masses is to change the kinematics; it has been shown that effects such as the 
apparent breaking of chiral symmetry can be neglected to good enough approximation for 
our purposes [53-55]. The thermal masses are given by: [55, 56] 



(B.l) 


(B.2) 


(B.3) 


where g 2 is the coupling constant of SU{2)l and gy is the coupling constant of U{1)y in 
the SM. We have neglected all SM Yukawa couplings except to the top quark, as well as 
contributions from A 3 . In our solutions we use = O.TIT, = 0.59T and rriL = 0.19T. 
In general, the asymmetry formed increases for large thermal masses. 

Interestingly, as the Y do not acquire a thermal mass, at temperatures above 1.4my^ 
the SM Higgs can decay into Y rather than the other way round (depending on the thermal 
mass of (p). While there could be concerns about the CP violation of these decays, these are 
ineffective as Hi is in equilibrium. This is in agreement both with ansatz calculations of the 
CP violation and the results from [54]. 

C Chemical Potentials 

Since (p, as well as the leptons, will be kept close to equilibrium above the EWPT it is 
unnecessary to have a Boltzmann equation for each particle species. Rather, we can solve 
our Boltzmann equations with the chemical potentials of these species, using (3.2) to obtain 
their number densities. This is where the advantage of using Maxwell-Boltzmann statistics 
lies: we can separate out the chemical potentials. Here we are not considering processes (or 
the relevant chemical potentials) below the EWPT. Our task is made simpler by the fact that 
all our chemical potentials can be written in terms of As our model is similar to standard 
treatments of B and L violation, such as leptogenesis, we can borrow the chemical potentials 
from [57], and just note that the segregation of B — L into cp is the only source of B — L 
violation. At the temperatures we will be considering, the population of H 2 is negligible and 
so does not affect the chemical potentials. By making the replacement 


B - L = 


(C.l) 
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we obtain the chemical potentials: 


(C.2) 

(C.3) 


-7 
4 

These can be related to the asymmetry generated in (p by using 

T2 

= -jH- (C^-4) 

To write this all in terms of the baryon asymmetry, we use 

28 

B = (C.5) 

D Extensions of the Neutron Portal 
D.l 2^3 scatterings 

We have considered the simplest gauge invariant combination oi B — L carrying SM particles, 
LH, but what about the second simplest? The next lowest dimensional gauge singlet B — L 
combination of SM particles is udd, the neutron portal. The neutron portal’s main issue 
is stability - the dark matter tends to decay. There are two ways to stabilise the neutron 
portal: temperature dependent couplings, as suggested in [10], or a symmetry. Both of 
these approaches require additional particle degrees of freedom. If we wish to impose a new 
symmetry, we must move beyond 2 —)• 2 scatterings. The essential neutron portal operator is 

XuRdp^dn, (D.l) 

where X is Dirac. It also possible to use left handed quarks, or for X to be Lorentz contracted 
with down type quarks, but there are no differences relevant to this discussion. 

Following the method of Section 2, the simplest way to make a stable version of the 
neutron portal is to introduce the effective operator 

gXuRdj^dna, (D.2) 

where X is a Majorana fermion and a is a complex SM singlet scalar.® As before, to satisfy 
(3.7) we require two generations of X. We will only consider one generation of quarks: as 
they remain in equilibrium any CP violation from having multiple generations will not change 
our discussion below. 

There are multiple ways to open up this effective operator, shown in fig. 9, but the 
only simple UV completion capable of generating an asymmetry is the last. In all others, 
there are rapid 2 —)• 2 flavour changing scatterings (fig. 10), which are problematic. As the 
2 —>■ 2 scatterings are mediated by only one heavy intermediate scalar, they dominate over the 
2 —>■ 3 scatterings. As CP violating effects involving the 2 —)> 3 scatterings must be as large 
as possible during freeze-out the delay of thermal non-equilibrium due to 2 —)• 2 scatterings 
results in a serious suppression of the final asymmetry. 
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Figure 9. Tree level UV completions of the stable neutron portal operator (D.2). 



Figure 10. Flavour changing scatterings. 


If the collision terms for u and the X are given by 

Wiu + Xi^d + d + a) = (l + ei)lFi, 

Wiu + X2^d + d + a) = {1 + €2)W2, (D.3) 

W{U + Xi^U + X2) = {l + €3)W3, 

then the relevant nnitarity constraints are given by 

eiWi + e2VF2 = 0, 

eiVFi + = 0, (D-4) 

62 ^ 1^2 ~ eahfs = 0. 

From this, it is clear that for any of these process to violate CP, flavonr changing scatterings 
mnst exist. Fnrther, to get signihcant CP violation withont delaying the departnre from 
thermal eqnilibrinm these scatterings mnst be comparable in size to the 2 —)• 3 scatterings. 
This argnment holds for any of the particles X conld scatter with; at least one flavonr chang¬ 
ing process mnst exist. As in this last UV completion there are no tree level scatterings to 
accomplish this, another process mnst be added.® 


®An alternative ADM scenario, using the neutron portal and CP violating decays, can be found in [58]. 
^Mediators for this process are contained in the other UV completions. 
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It seems that the only way to obtain a working ADM nentron portal using 2 —>■ 3 scatter¬ 
ings is to have correlated couplings, which is not obviously an improvement over temperature 
dependent couplings. Any full theory of 2 —)• 3 scatterings must explain why two seemingly 
unrelated processes conspire to allow an asymmetry to be formed. While the neutron portal is 
the second simplest SM gauge singlet operator involving baryon or lepton number, it actually 
requires a complex description. We must look elsewhere for a generic method for ADM via 
scatterings. 

D.2 What about the symmetric component? 

As discussed in Section 2.1, no existing models can use the same process to annihilate the 
symmetric component and generate an asymmetry but is it possible to design one that does? 
Any model satisfying this condition must have four properties: 

1. The particles going out-of-equilibrium must be stable. 

2. The process must freeze-out late enough for the symmetric component to be annihilated, 

?freeze-out < mDM/25. 

3. Dark matter must couple to quarks, as a lepton asymmetry cannot be reprocessed below 
Tewpt and ADM requires that moM is below the EW scale. 

4. CP violation must be significant at freeze-out, at least 10“® (a conservative lower bound 
by analogy with [18], which used kinematic suppression of the washout to enhance the 
asymmetry formed). 

As there are strict bounds on light non-SM coloured particles [59], DM must couple to a gauge 
singlet combination of quarks carrying baryon number, so some extension of the neutron por¬ 
tal is necessary. At least one of the mediators for any interaction involving the neutron portal 
must be coloured, making it potentially detectable at the LHC. We can derive an upper 
bound on m-mediators from properties 2 and 4. As dark matter is very light (order GeV), this is 
in tension with collider searches, which require that new coloured particles are heavier than 
about a TeV (for gluino-like particles) and 440 GeV (for single light flavoured squark-like 
particles) [59]. Models with small mass splittings can weaken these bounds, however as we 
are considering ADM the dark matter will always be significantly lighter than the coloured 
mediators. 

The simplest model which satisfies these conditions is (D.2),^’^ so we can get the most 
forgiving bounds from it. More complicated extensions of the neutron portal will in general 
have tighter bounds, as additional mediators suppress the interaction. To simplify matters 
we will consider all mediators to have the same mass; light SM gauge singlet mediators tend 
to lead to rapid D conserving scatterings amongst the dark sector particles. We estimate 
the freeze-out temperature by comparing the interaction rate T to the Hubble time: freeze- 
out occurs when T ~ Lf. These results agree with the numerical calculations of the similar 
neutrino portal model. Requiring that freeze-out occurs late enough (point 2 above) gives 

^mediator ^ 150 X A GeV, (D-S) 

we make X Dirac, rather than Majorana, and kinematically disallow decays dark matter will be an 
admixture of X and cr, similar to [58]. 
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where A is the couplings strength, assuming all couplings are the same. The CP violation at 

r^3 

low temperatures is bounded by a loop factor ——. From 4 we get 

^^mediator 

^TTmediator < 200 X A GeV. (D. 6 ) 

While limits on coloured particles are model dependent, coloured particles of order 150 GeV 
would almost certainly have been seen at the LHC. The squark limits can be avoided for 
A > 3, but this raises the issue of perturbativity and low energy Landau poles. So strong a 
Yukawa interaction between quarks and the mediating particles means that there would be 
a very strong detection possibility via, for example, monojet searches. This is an interesting 
avenue for future work. 

Unfortunately it seems that no operator allows one to use a single scattering process to 
annihilate the symmetric component and generate an asymmetry through asymmetric freeze- 
out without using couplings greater than one. Further, this requires at least 2 —)• 3 scatterings, 
so for the same process to annihilate the symmetric component and generate an asymmetry 
simultaneously, it must rely on a conspiracy with flavour violating scatterings. 


E Other 2 —)■ 2 ADM Scattering Operators 

As was shown in Appendix D, 2 —)• 2 scatterings have signihcant advantages over 2 —)• 3 
scatterings. Other than the neutrino portal, what other 2 — 7-2 scatterings are possible? If 
the SM particles involved are not required to form a gauge singlet then there are dozens of 
operators that can be used. Due to the sheer number of operators (many being only trivially 
different) an exhaustive list would be exhausting both to the reader and the authors. Instead, 
we will discuss the possibilities using an example. 

For the example, consider WIMPy baryogenesis [17]. WIMPy baryogenesis also uses 
GP violating scatterings to generate a baryon asymmetry, albeit using the "WIMP miracle" 
instead of an asymmetry to set the relic abundance of dark matter. Despite this, it is possible 
to tweak WIMPy baryogenesis so that it becomes a theory of ADM. Fortuitously, the different 
ways of modifying WIMPy baryogenesis reveal all the salient model building concerns of 
2 —)• 2 ADM scattering operators. For stock WIMPy baryogenesis, the relevant additions to 
the Lagrangian are [17] 

AT = gL'il^X'^X -h XHi^f + H.c, (E. 1) 

where X and / are singlet Dirac fermions, the dark matter, and ^ carries the same quantum 
numbers as L. In the original model, the X freeze-out, not carrying any baryon number, 
behave as WIMPs. GP violation in these scatterings stores lepton number in the -0. The 
ips (which are heavier than X) decay via ^|J — H*f. Dangerous decays giving lepton num¬ 
ber back to the visible sector are forbidden by a Z 4 symmetry. Of course, there are other 
WIMPy baryogenesis models with, for example, X coupling to quarks but these give a similar 
story. There are two different approaches to turning the WIMPy baryogenesis operator into a 
model of ADM, corresponding to the two broad categories of 2 —?■ 2 ADM scattering operators. 
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The first approach corresponds to a class of theories similar to the toy model of [12].^^ 
In this first class of theories, particles that store an asymmetry go out-of-equilibrium and then 
decay into lighter particles, transferring the asymmetry. This decay is required to allow the 
annihilation of the symmetric component. To accomplish this, and to get a tighter relationship 
between dark and visible matter, we can do away with X. There is already a conserved U(l) 
symmetry between L and il), so the process 

LL (E.2) 

could create an asymmetry through asymmetric freeze-out, with ■0 going out of equilibrium. 
The asymmetry stored in ^ is transferred to / as in WIMPy baryogenesis, but now dark mat¬ 
ter consists solely of /. In place of X, a mechanism to annihilate the symmetric component 
of / is necessary but for essentially the same degrees of freedom we now have Odm — 
naturally. 

In the second approach, which is similar to the neutrino portal model and [13], particles 
that do not carry B — L — D depart from equilibrium. In the WIMPy baryogenesis scenario, 
we can replace one of the Xs in this operator with a Majorana fermion Y, giving 

AC = gL^pX^Y + XHPpf + H.c. (E.3) 

Now X is free to be light, and to remain in thermal equilibrium. As in the neutrino portal 
model, the Majorana fermion Y departs from equilibrium, and subsequently decays. Unlike 
WIMPy baryogenesis, we insist that X carry a baryon number equivalent. This leads to 
two component dark matter, consisting of X and /. In order to ensure the stability of / we 
introduce a Z 2 symmetry, with Y, ip and / being the negative parity states. If X is lighter 
then /, Nx — (Nf + N^) conservation ensures the stability of X. To annihilate the symmetric 
component, it is possible to gauge U(l)jv_, 5 .-(Arj.+Ar,^)- This should produce a similar result to 
the neutrino portal as the Boltzmann equations and types of CP violation are quite similar, 
though it is non-minimal as it only involves one SM particle. 

All other operators using this mechanism of asymmetric freeze-out should fall broadly 
into one of these two categories. While the neutrino portal is the minimal case, adding 
complexity to the dark sector can lead to some potentially interesting phenomenology, such 
as dark forces. 
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